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Abstract 

We propose two ansatze of lepton mass matrices with four texture zeros, and confront 
them with current experimental data on neutrino oscillations. The parameter space of 
each ansatz is carefully explored. We find that both ansatze can accommodate the normal 
hierarchy of neutrino masses and the bi-large pattern of lepton flavor mixing. Their 
predictions for the effective mass of the tritium beta decay and that of the neutrinoless 
double beta decay are too small to be detectable, but leptonic CP violation at the 
percent level is allowed. Some discussions are also given about the seesaw invariance of 
the four-zero texture of Dirac and Majorana neutrino mass matrices. 
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1 Introduction 



The recent KamLAND [1] and SNO [2] experiments have provided us with very compelling 
evidence that the solar neutrino deficit is due to the matter-enhanced v e — > oscillation with 
a large mixing angle 6 snn ~ 32°. Meanwhile, the K2K [3] and Super-Kamiokande [4] experi- 
ments have convinced us that the atmospheric neutrino anomaly is attributed to the — > v T 
oscillation with another large mixing angle # atm ~ 45°. In contrast, the non-observation of the 
V e — > V e oscillation in the CHOOZ experiment [5] indicates a rather small (even vanishing) 
mixing angle 9 chz < 13°. These three mixing angles may straightforwardly be related to the 
elements of the 3x3 lepton flavor mixing matrix V, which links the neutrino mass eigenstates 
(ui, v-ii ^3) to the neutrino flavor eigenstates (z/ e , z/ M , v T \ To a good degree of accuracy, 6 snn and 
#atm describe the flavor mixing effect between the 1st and 2nd lepton families and that between 
the 2nd and 3rd lepton families, respectively; while the small mixing angle # c h z is responsible 
for the flavor mixing effect between the 1st and 3rd lepton families. Thus the lepton flavor 
mixing matrix V performs a bi-large mixing pattern, which is quite different from the tri-small 
mixing pattern of the quark flavor mixing matrix [6]. 

To interpret the bi-large lepton flavor mixing pattern, many phenomenological ansatze of 
lepton mass matrices have recently been proposed [7]. A particularly interesting category of 
the ansatze focus on texture zeros of charged lepton and neutrino mass matrices in a given 
flavor basis from which some nontrivial relations between flavor mixing angles and lepton 
mass ratios can be derived. The typical example is the Fritzsch ansatz [9] for symmetric lepton 
mass matrices, 

/ Ci, v \ 
M hv = Q, v B ltV , (1) 
V B liV A^j 

in which six texture zeros are included §. It has been shown by one of the authors [11] that 
this type of lepton mass matrices can naturally predict a normal hierarchy of neutrino masses 
and a bi-large pattern of lepton flavor mixing angles. Furthermore, Fukugita, Tanimoto and 
Yanagida [12] have demonstrated that very similar phenomenological predictions can also be 
achieved from a simple but interesting ansatz of lepton mass matrices based on both the 
Fritzsch texture and the seesaw mechanism [13]. 

The present paper aims to analyze the generalized Fritzsch ansatz of lepton mass matrices 
with four texture zeros, 




M ltV = C hv B hv B ltV , (2) 



and its consequences on the neutrino mass spectrum, flavor mixing and CP violation. It is 
well known that the four-zero texture of quark mass matrices is more successful than the 
six-zero texture of quark mass matrices to interpret the strong hierarchy of quark masses and 
the smallness of flavor mixing angles. The spirit of lepton-quark similarity motivates us to 
conjecture that the lepton mass matrices might have the same texture zeros as the quark mass 
matrices. Such a conjecture is indeed reasonable in some specific models of grand unified 



■l- For instance, a lot of interest has recently been paid to possible two-zero textures of the neutrino mass 
matrix in the flavor basis where the charged lepton mass matrix is diagonal [8] . 

§A pair of off-diagonal texture zeros of the charged lepton (Mi) or neutrino (M v ) mass matrix have been 
counted, due to symmetry, as one zero [10]. 
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theories [7], in which the mass matrices of leptons and quarks are related to each other by a 
new kind of flavor symmetry. That is why the four-zero texture of lepton mass matrices has 
been considered as a typical example in some model-building works [14]. However, a careful 
and complete analysis of its phenomenological implications has not been done in the literature. 

Naively, there is no doubt that the four- zero texture of M\ %v in Eq. (2), which has more 
free parameters than the Fritzsch texture of in Eq. (1), must be able to interpret the 
observed bi-large pattern of lepton flavor mixing. To improve the analytical calculability and 
numerical predictability, one may follow a realistic strategy to concentrate on part of the 
whole parameter space of the four-zero texture of lepton mass matrices. The same strategy 
has actually been adopted in the study of the four-zero texture of quark mass matrices [15]. 
In this paper, we shall consider two simplified versions of Mi and M v given in Eq. (2): 

1. \Bi\ — m M and \B U \ = m 2 , where m M and m 2 stand respectively for the physical masses 
of /j, and u 2 . This interesting case, to be referred to as Ansatz (A), has been briefly 
discussed in Ref. [14]. 

2. \B t \ = \B t \ and \B U \ = \A U \. This specific case, to be referred to as Ansatz (B), has not 
been discussed in the literature. 

Our main purpose is to explore the allowed parameter space of each ansatz and its implications 
on the neutrino mass spectrum and lepton flavor mixing measurables. We find that both 
ansatze are favored by current experimental data on neutrino oscillations. Their predictions 
for the effective mass of the tritium beta decay and that of the neutrinoless double beta decay 
are too small to be detectable, but leptonic CP violation at the percent level is definitely 
allowed. Finally, we present some brief discussions about the seesaw invariance of the four- 
zero texture of Dirac and Majorana neutrino mass matrices. 



2 Framework 

Note that all non-zero elements of Mi and M u in Eq. (2) are in general complex. If the 
condition 

arg(A^) + axg(B ltV ) = 2 axg(S, )I/ ) (3) 
is satisfied, then M^ v can be decomposed as 

Mi = PjMiPi, M U = P?M V P U , (4) 

where 



Pl,v — 






\Cl,v\ 







\Bi,v\ 







\Bi,v\ 




id, 











iB, 

e l ' v 














(5) 



with arg(A^) = 2 7z iy , arg(£^) = 7^, arg(i^) = 2/^„ and arg(C til/ ) = a l>v + P l>v . 

The real symmetric matrices Mi and M v can be diagonalized by use of the following unitary 
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transformations: 
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( m l 
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XJ T V M V U V = 





m 2 







1 o 





m 3 / 



(6) 



in which (m e ,m^,m T ) and (mi,m2,m 3 ) denote the physical masses of charged leptons and 
neutrinos, respectively. The lepton flavor mixing matrix V arises from the mismatch between 
the diagonalization of the charged lepton mass matrix Mi and that of the neutrino mass matrix 
M v . Therefore, we obtain V = Uf (P*P U ) U*, whose nine matrix elements read explicitly as 



y — TJl Tjv* p ia , Tjl ttv* , Tjl ttv* if 
Vpq — U lv U lQ e l - J 2p U 2o e U 3v U 3a e 5 



jl j jv* la 
'lp U lq ' 



(7) 



where the subscripts p and q run respectively over (e, //, r) and (1,2,3), and the parameters 
a, P and 7 are defined as 



a = a v - a t , 



P = - A , i = i v -ii 



(8) 



Note that the overall phase of V has nothing to do with the experimental observables. Hence 
only two combinations of three phases (a, (3, 7) are physically relevant. For simplicity, we take 
7 = in the following. 

The matrix elements of V depend both on the ratios of lepton masses, 



m M m T m 2 m 3 



(9) 



and on the phase parameters a and (5. As the values of m e , m M and m T have precisely been 
measured [16], we have x l w 0.00484 and y l ps 0.0594 to a good degree of accuracy. The other 
four free parameters can be determined or constrained from the present experimental data 
on neutrino oscillations. As the neutrino mass-squared differences of solar and atmospheric 
neutrino oscillations are given by 



(10) 



Am L = 


2 2 

m 2 — m 1 


to to 


1 - 


x v 


Am atm = 


2 2 
m 3 — m 2 


= m l 


1 - 


-vl 



the observed hierarchy Am 2 un Am^ tm may impose a very strong constraint on the values 
of (x v ,y v ): 



Am 2 



vl 



atm 



X, 



yl 



< 1 



On the other hand, the mixing factors of solar, atmospheric and CHOOZ reactor neutrino 
oscillations are related to the matrix elements of V in the following way: 



sin 2 26 s 
sin 2 1Q R 

sin 2 26 chz 



' atm 



4|v; 1 | 2 |\/ e2 

4|V^ 3 | 2 (l 

4|K 3 | 2 (1 - |v; 3 | 2 ) 



;i2) 



In view of the KamLAND [1], SNO [2], K2K [3], Super-Kamiokande [4] and CHOOZ [5] data 
on neutrino oscillations, we have Am s 2 un G [5.9, 8.8] x 10~ 5 eV 2 , sin 2 6» sun G [0.25, 0.40] [17]; 
Am 2 tm G [1.65, 3.25] x 10~ 3 eV 2 , sin 2 20 atm G [0.88, 1.00] [18]; and sin 2 26 chz < 0.2 at the 90% 
confidence level. With the help of these experimental results, the allowed ranges of x„, y„, a 
and P can be obtained from Eqs. (11) and (12). 

Once the values of x v and y v are determined or constrained from current experimental 
data, we are able to calculate the absolute values of three neutrino masses by use of Eq. (10): 



m 3 



m 2 



mi = 



1 


>- 


vl\ 


Vv 




>- 

1 


vl\ 


>- 


x 2 \ 


x v 





Am 2 



atm ) 



Ami 



v^ 1 



Am sun 



(13) 



Xf 



In addition, interesting predictions can be achieved for the effective mass of the tritium beta 
decay (m) e and that of the neutrinoless double beta decay (m) ee : 



HI EE EHl^l 2 ) 
i=l 

(xly 2 u \V el \ 2 + yl\V e2 \ 2 + \V e3 \ 2 ) , 

E (^) 



(m) 



= m 3 \x v y v V 2 x + y v V 2 2 + 



(14) 



The present experimental upper bound on (m) e is (m) e < 3 eV [16], while the sensitivity of the 
proposed KATRIN experiment is expected to reach (m) e ~ 0.3 eV [19]. In comparison, the 
upper limit (m) ee < 0.35 eV has been set by the Heidelberg-Moscow Collaboration [20] at the 
90% confidence level *. The sensitivity of the next-generation experiments for the neutrinoless 
double beta decay is possible to reach (m) ee ~ 10 meV to 50 meV [22]. 

The strength of CP violation in neutrino oscillations, which is measured by the Jarlskog 
invariant J [23], can also be predicted from the four-zero texture of lepton mass matrices 
under consideration. Indeed, J is defined through the following equation: 



c,k 



(15) 



where the subscripts (a,b,c) and (i,j,k) run respectively over (e, //, t) and (1,2,3). The 
magnitude of J depends both on (x u ,y u ) and on (a,/3). If \ J\ ~ 1% is achievable, then 
leptonic CP- and T-violating effects could be measured in a variety of long-baseline neutrino 
oscillation experiments [24] in the future. 

'if the reported evidence for the existence of the neutrinoless double beta decay [21] is taken into account, 
one has 0.05 eV < (m) ee < 0.84 eV at the 95% confidence level. 
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3 Ansatz (A) 



Now let us consider Ansatz (A), in which the requirements \Bi\ = and \B V \ = m 2 are 
imposed on Mi and M v in Eq. (2). Similar conditions (\B U \ = m c and \B&\ = m s ) have 
actually been taken in some literature [14] for the four-zero texture of quark mass matrices. 
Following Eqs. (3) and (4), we factor out the complex phases of M^ v . Then three elements of 
the real symmetric mass matrix Mi u can be expressed in terms of its three mass eigenvalues: 



\M 

W 

\c,\ 



m T — m e , 



m e m T {m T — m e — to m ) 

m T — m e 
m e m^m T \ 1 / 2 



1/2 



(16) 



and 



1713-171! , 



\B V 



m 1 m 3 ( y m 3 — m\ — m 2 ) 



/ m 1 m 2 m 3 \ 
\m 3 — mi) 



m 3 — mi 

1/2 



1/2 



(17) 



The elements of the unitary transformation matrix Ui yU , which is used to diagonalize M^ u in 
Eq. (6), can in turn be expressed in terms of the ratios x l v and y l u as follows (the indices "Z" 
and 'V are neglected for simplicity): 



U n 
U 12 
U 13 
U 2 i 
U22 

u 23 
U31 
u 32 
u 33 



+i 



= + 



1 



1/2 



(l + x)(l-x 2 y 2 ) 

x(l-y- xy) 

[1 + x)(l - y)(l - xy) 
1/2 



1/2 



+ 



x 2 y 3 



I — y)(l — x 2 y 2 
x 

(l + x)(l + xy) 



1/2 



+ 



= + 



1/2 



1 1/2 



+i 



1/2 



+ 



l-y-xy 
(l + x)(l-y) 

xy 

(l-y)(l + xy) 

x 2 y(l -y-xy) 
(1 + x)(l - x 2 y 2 ) 

xy 

(l + x)(l-y)(l-xy) 

1 — y — xy 

(l-y)(l-xV) 



1/2 



(18) 
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Note that Un (for % = 1,2,3) are imaginary, and their nontrivial phases are due to the negative 
determinant of M^ v . 

The four free parameters x v , y u , a and (5 can be constrained by use of Eqs. (11) and (12) 
as well as current data on neutrino oscillations. Their allowed ranges are shown in Fig. 1. We 
see that x v ~ 0.86 and y v ~ 0.35 typically hold. Thus the neutrino mass spectrum satisfies 
nil < m 2 < m 3- For a and j3 varying from to 2tt, we find that about half of the whole (a, (3) 
parameter space can be excluded. Note that the correlation between a and (5 in Ansatz (A) is 
not as strong as the phase correlation in the Fritzsch ansatz of lepton mass matrices [11, 12]. 
The reason is simply that the contribution of Mi to V is much smaller in Ansatz (A) than 
in the Fritzsch ansatz. Hence the relative phases between Mi and M v in the former cannot 
significantly affect the magnitudes of nine matrix elements of V. 

Fig. 1 also shows the ouputs of sin 2 2# atm versus sin 2 9 sun and sin 2 2# c h z versus R u restricted 
by Ansatz (A). It can be seen that larger values of sin 2 9 snn roughly correspond to smaller 
values of sin 2 2# atm . In addition, the ansatz predicts sin 2 26> c h z > 0.08, a lower bound which 
is easily accessible in the upcoming long-baseline neutrino oscillation experiments [24]. If the 
upper limit sin 2 29 chz < 0.1 instead of sin 2 29 chz < 0.2 is input in the numerical calculations, 
one will arrive at sin 2 2# atm < 0.91. Such a low value of sin 2 2# atm is tolerable, but not favored 
by current data. It becomes clear that the mixing angles 9 sun , 9 atm and # c h z are strongly 
correlated with one another in Ansatz (A). Thus more precise data on three mixing angles 
may provide a sensitive test of this phenomenological scenario. 



The result y 2 v ~ 0.1 implies that m 3 ps yAm atm is an acceptable approximation. More 
exactly, we obtain m 3 ps (4.3 - 6.1) x 10 -2 eV, m 2 ~ (1.4 - 2.3) x 10~ 2 eV and m x 
(1.1 — 2.1) x 10~ 2 eV from Eq. (13). In calculating the allowed ranges of m\ and m 2 , we have 
ignored their correlation induced by the model itself. This generous treatment has no conflict 
with the plot of (x u , y u ) in Fig. 1, in which x v < 1 results from the correlation between m 1 /m2 
and m 2 /m 3 . The sum of three neutrino masses is consistent with m\ +m 2 + m 3 < 0.71 eV, an 
upper bound set by the recent WMAP data [25]. We compute the effective mass of the tritium 
beta decay and that of the neutrinoless double beta decay by use of Eq. (14), and present the 
numerical results in Fig. 1. One can see that (m) e ~ 10~ 2 eV and (m) ee ~ 10~ 3 eV typically 
hold. Both quantities are too small to be measured in practice. Similarly, there is no hope to 
detect the effective (kinematic) masses of muon and tau neutrinos [26]. The numerical results 
for the Jarlskog parameter J and the smallest matrix element |T4 3 | are also shown in Fig. 1. 
We find that the magnitude of J may nearly be 1%, if |V e3 | is close to its upper bound. It is 
possible to measure leptonic CP violation of this order in the future neutrino factories, if the 
terrestrial matter effects can be under control. 

Finally we illustrate the typical texture of M^ u by taking x u ~ 0.86, y u ~ 0.35 and 
m 3 ps 0.05 eV. The result is 



Mi ~ 1.78 GeV x 



M u ~ 3.50 x 10~ 2 eV x 



/ 0.0041 

0.0041 0.059 0.016 

V 0.016 1 

/ 0.56 \ 

0.56 0.50 0.55 

V 0.55 1 / 



(19) 



It becomes obvious that lepton flavor mixing is dominated by the neutrino sector, as the 
matrix elements of Mi have a very strong hierarchy. 
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4 Ansatz (B) 



We proceed to consider Ansatz (B), in which the requirements \Bi\ — \Bi\ and \B V \ = \A V \ 
are imposed on Mi and M v in Eq. (2). Note that the condition \B t \ = \Bi\ is similar to 
\B U \ « \B U \ and \B^\ ~ \B^\ for the four-zero texture of quark mass matrices [15], in view of 
the fact that charged leptons have a strong mass hierarchy as quarks. Because the condition 



Bi\ leads to \B. 



\Bi\ 

to the condition IS 



'i\ in the leading-order approximation, it is essentially equivalent 

'i\ = taken in Ansatz (A). In contrast, the requirement \B V \ = \A U \ is 



motivated by the experimental fact that the mixing angle of atmospheric neutrino oscillations 
is about 45° [4] (namely, the diagonalization of the (2,3) subsector of M v may give rise to 
a rotation angle of 45°, if the condition \B V \ = \A U \ is satisfied). Such a phenomenological 
hypothesis for the texture of M v results in an apparent "structural asymmetry" between M\ 
and M„, but its consequences on lepton flavor mixing are simple and interesting. Following 
Eqs. (3) and (4), we can factor out the complex phases of M^ v . Although we are able to 
exactly express both (\Ai\, \Bi\, \Ci\) and (\A U \, \B V \, \C V \) in terms of the corresponding mass 
eigenvalues, the formulas for the former are too complicated to be instructive. It is therefore 
better to make some analytical approximations in deriving \A{\, \B{\ and \C{\. In view of 
the strong mass hierarchy in the charged lepton sector, we expect that \A t \ ^> \Bi\ ^> \C{\ 
naturally holds. Then we obtain 



\M 



\Bi\ 



\Ci\ 



m, 




(20) 



to a good degree of accuracy. In contrast, the expressions of 

m 3 + m 2 -mi 



\B U \ and |CJ are exact: 



\BJ = - 



(777,3 + m 2 + mi) (m 3 — m 2 — mi) (7773 — m 2 + mi) 



m-i + m 2 — mi 



1/2 



|C„| 



2mi777 2 7773 \ 



(21) 

,7773 + 7772 — mi) 

The unitary transformation matrix Ui, which has been defined to diagonalize M; in Eq. (6), 
can approximately be given as 



U, 



x, 



IX, 



Obi 



fx, 



V w 



1 _ ZL - £L 



-Vl 



Vl 



Vl_ 
2 



(22) 



where the fact of x l ~ yf has been taken into account. In addition, the matrix elements of U v 
can be expressed in terms of x u and y v as follows: 

nl/2 



+i 



(1 + y v + x v y v ) 



(1 + x u )(l + x v y v )(l +y v - x v y v ) 
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tf 

u 12 



TF 

u 13 



TF 

u 21 



TF 

u 22 



TF 

u 23 



TF 

u 31 



TF 

u 32 



TF 

u 33 



+ 



= + 



—I 



+ 



= + 



= +i 



= + 



•^1/(1 Vv -^uVu) 



,1/2 



(1 + x v )(\ - y v ){\ + y v - x v y v ) 

XyVlil -yv + x u y v ) 
(1 - y v )(l + x v y v )(l + y„- x v y v ) 

1 1/2 

x v (l + y v + x u y v ) 
2(l + x u )(l + 

(l-y„- x v y v ) 



1/2 



nl/2 



2{l + x v ){l-y v )_ 

(l-y u + x v y v ) 
2(1 -y v ){l + x v y v ) 



x v {\ -y v - x u y u )(l -y u + x v y v ) 
2(1 + x v )(l + x u y v )(l +y„- x v y v ) 

(1 + y u + x„y v )(l -y„ + x u y u ) 



1/2 



1/2 



2(1 + x v ) (1 - y v ) (1 + y v - x v y v ) 

(l+y v + x v y v ){\ -y u - x v y v ) 
2(1 - y v ){l + x v y v )(l + 

Vv x v y v ) 



1/2 



(23) 



Note again that U l a and U\ x (for i — 1, 2, 3) are imaginary, and their nontrivial phases are due 
to Det(M,) < and Det(M^) < 0. 

With the help of Eqs. (11) and (12), one may compute the ranges of x v , y u , a and (3 allowed 
by current data on neutrino oscillations. The numerical results are shown in Fig. 2. We see 
that the allowed range of x v in Ansatz (B) is much larger than that in Ansatz (A). As both 
x v < 1 and y v < 1 hold, the neutrino mass spectrum satisfies mi < m 2 < m 3 . Note that there 
is little restriction on the phase parameters at and (5. This feature can easily be understood: 
the contribution of M v to V is dominant over that of M t to V, thus the magnitudes of nine 
matrix elements of V are essentially insensitive to the relative phases between Mi and M v . 

The outputs of sin 2 2# atm versus sin 2 # sun and sin 2 20 c h z versus R v are also shown in Fig. 2. 
We see that sin 2 26* atm > 0.94 holds and sin 2 2# atm 1 is particularly favored. The latter is a 
natural consequence of the specific texture of M u , whose (2,3) subsector can be diagonalized 
by a rotation of 45°. Thus |V^ 3 | l/V^ leads to sin 2 29 atm m 1 in Ansatz (B). Another 
feature of Ansatz (B) is that changing the upper bound of sin 2 26 chz from 0.2 to 0.1 does not 
significantly affect the allowed range of sin 2 26 l atm . However, the correlation between sin 2 29 chz 
and R v is stronger in Ansatz (B) than in Ansatz (A). 

The approximation m 3 ps \J Am atm is reasonably good in Ansatz (B). Numerically, we 
obtain m 3 ^ (4.2 - 5.8) x 10~ 2 eV, m 2 w (0.8 - 1.4) x 10" 2 eV and m x ps (0.3 - 1.1) x 10~ 2 eV 
from Eq. (13). This neutrino mass spectrum is quite similar to that in Ansatz (A). The results 
for the effective mass of the tritium beta decay and that of the neutrinoless double beta decay 
are also illustrated in Fig. 2. One can see that (m) e ~ 10~ 2 eV and (m) ee ~ 10~ 3 eV hold in 
Ansatz (B). Both of them are too small to be detected in reality. In addition, Fig. 2 shows 
that the magnitude of the Jarlskog parameter J may nearly be 1.5%, if jV^] is larger than 
0.15. This result implies that it is possible to observe leptonic CP violation in the future 
long-baseline neutrino oscillation experiments. 
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To illustrate the texture of Mi tU , we typically take x v ~ 0.6, y v ~ 0.21 and m 3 ss 0.05 eV. 
The numerical result is 



Mi 



My 



0.0042 \ 

1.77 Go\ x j 0.0042 0.064 0.064 

0.064 1 / 

/ 0.41 



2.71 x 10 _i eV x 



0.41 1 0.80 
V 0.80 1 



(24) 



The similarity and difference between two ansatze are therefore obvious. 



5 Seesaw 



Two ansatze of lepton mass matrices discussed above can self-consistently describe the ob- 
served features of lepton flavor mixing, but they give no interpretation about why the masses 
of three neutrinos are so tiny. A simple way to improve our phenomenological ansatze is to 
incorporate them with the elegant idea of seesaw [13]. In the seesaw mechanism, the small- 
ness of left-handed Majorana neutrinos is attributed to the existence of heavy right-handed 
Majorana neutrinos, 

M v « M^M^Ml , (25) 

where Md and Mr denote the Dirac neutrino mass matrix and the heavy Majorana neutrino 
mass matrix, respectively. In some grand unified theories (such as the SO (10) model [27]), one 
takes [Md,M u ] = 0, where M u represents the up- type quark mass matrix. The mass matrix 
Mr is practically unknown in almost all reasonable extensions of the standard model. Hence 
specific textures of M R and M D have to be assumed, in order to determine the texture of M v . 
Given M u and M D , on the other hand, one can calculate Mr by use of Eq. (25): 



M 



M^M^M 



D 



(26) 



(27) 



The scale of Mr stands for the scale of new physics in this simple seesaw picture. 

To be specific, we assume that M D = M u holds and it has the same texture zeros as M 
and My have: 

/ o C u 

M D = C u B u B u 
V B u A U: 

Only if the condition DetM, ^ is guaranteed for My in Eq. (2), one may obtain the inverse 
matrix of My as follows: 

/ By — AyBy AyCy ~ B yC y \ 



o el 



M7 1 = 



\ 



AyCy 
-ByCy 



(28) 



Then the texture of M R can be determined from Eq. (26) with the help of Eqs. (27) and (28): 



(29) 





( 


C R 


\ 


M R = 




Br 


Br 




V o 


Br 


Ar] 
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where 



A 2 

A _ U 

— ~r > 



A 



V 



^ A n B u B U C U A u C n B u 

A v C u A V C V 

~ 2B n C u C\B V 2B U C U B V C\B 2 V 

R ~~ 4 r r* 2 at ac 2 ' 1 

l\ v Ly u Kj v 1\ v Kj v 1\ v Vj v 

Cn = (30) 

We see that the texture zeros of Md and M v manifest themselves again in Mr, as a consequence 
of the inverted seesaw relation given in Eq. (26). Therefore, all four lepton mass matrices 
(M, Md, Mr and M v ) are structurally parallel to one another. Such a structural similarity of 
lepton mass matrices, which is seesaw-invariant, might follow from a universal flavor symmetry 
hidden in a more fundamental theory of fermion mass generation. In particular, the underlying 
flavor symmetry must be related to the texture zeros of lepton mass matrices. It is worth 
mentioning that the Fritzsch texture of lepton mass matrices in Eq. (1) does not have the 
interesting property of seesaw invariance. Thus we argue that the four-zero texture of lepton 
mass matrix might be more attractive for model building at the energy scale where the seesaw 
mechanism works. 

Now let us give an order-of-magnitude estimate of the matrix elements of M R by taking 
the following phenomenologically-favored pattern of M u [15]: 

/ ^m u m c \ 



M u ~ 



(31) 



\ y/m u m t m t J 

in which the relevant complex phases have been neglected for simplicity. Typically taking 
m u /m c ~ m c /m t ~ 0.0031 and m t m 175 GeV at the electroweak scale [16], we obtain 

/ 5.4 x 10~ 8 

M R ~ 8.8 x 10 14 GeV x 5.4 x 10~ 8 9.6 x 10~ 6 3.1 x 10~ 3 | (32) 

V 3.1 x lO" 3 1 

from Eq. (19) for Ansatz (A); and 

7.3 x 10- 8 

Mr ~ 1.1 x 10 15 GeV x I 7.3 x 10~ 8 9.6 x lO" 6 3.1 x 10~ 3 I (33) 





3.1 x 10~ 3 

from Eq. (24) for Ansatz (B). We see that the structure of Mr is strongly hierarchical in 
either case. The scale of M R is about 10 15 GeV, close to the typical scale of grand unified 
theories A GU t ~ 10 16 GeV. 

It is worth remarking that the phase parameters of M D and M u have been ignored in 
estimating the matrix elements of Mr. If those complex phases are included, it is possible to 
get CP violation in the lepton-number- violating decays of heavy Majorana neutrinos [28]. 
However, a successful interpretation of the observed matter-antimatter asymmetry of the 
universe via the leptogenesis mechanism [28] is rather nontrivial, because the details of M D 
and Mr have to be taken into account. Further discussions on this topic are interesting but 
beyond the scope of this paper. 
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6 Summary 



We have proposed and discussed two phenomenological ansatze of lepton mass matrices with 
four texture zeros. The parameter space of each ansatz has been carefully analyzed by use of 
current experimental data on neutrino oscillations. We demonstrate that the normal hierarchy 
of neutrino masses and the bi-large pattern of lepton flavor mixing can be accommodated in 
both ansatze. Their predictions for the effective mass of the tritium beta decay and that of 
the neutrinoless double beta decay are too small to be detected in practice. However, we find 
that leptonic CP violation at the percent level is possible for either ansatz. The correlation of 
relevant observable quantities in each ansatz allows us to test its validity, once more accurate 
experimental data become available. This property may also allow us to distinguish between 
these two different ansatze. 

For the purpose of illustration, we have presented some brief discussions about the seesaw 
realization of our phenomenological scenarios. It is clear that the existence of heavy right- 
handed Majorana neutrinos at the scale of 10 15 GeV or so may naturally interpret the smallness 
of left-handed Majorana neutrino masses. This observation would be useful for model building, 
from which some deeper understanding of the neutrino mass generation and lepton flavor 
mixing could be gained. 

One of us (H.Z.) is grateful to the theory division of IHEP for financial support and 
hospitality in Beijing. This work was supported in part by the National Natural Science 
Foundation of China. 
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Figure 1: The parameter space and phenomenological predictions of Ansatz (A). 
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Figure 2: The parameter space and phenomenological predictions of Ansatz (B). 
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